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1. INTRODUCTION 

Due to advantages of superior massive connectivity and spectral efficiency, non-orthogonal multiple 
access (NOMA) has received considerable consideration as the prominent approach for future wireless net- 
works [1], [2]. As promising feature of NOMA, the same radio resources with different power factors based on 
their channel conditions share for signals from the plurality of users. At the receiving end, the user with poor 
channel conditions is able to decodes its own message as considering other user’s messages as interference. In 
constrast, the users with better channel conditions need successive interference cancellation (SIC) to aliviate 
another users’ messages and then they decode their own signal [3], [4]. By allocating different levels of power 
to users, NOMA relying on power domain to permit multiple users to achieve the same non-orthogonal re- 
sources [5]-[8]. The users’ messages would be detached with the help of SIC. In such, improving connectivity 
in IoT needs implementation of NOMA technique [9]. In [10] the authors targeted to security-required and the 
delay-sensitive users on the same non-orthogonal resource to satisfy requirements in IoT employing NOMA. 
They proved that the proposed system relying on NOMA paradigm is better than the benchmark OMA scheme. 
The throughput and energy efficiency are studied in the NOMA scenarios to perform a cellular massive IoT in 
[11]. further improving network connectivity, NOMA can be applied with cognitive radio and millimeter wave 
networks [12], [13]. In the maner of cooperative scheme, NOMA exhibits some advantages in term of outage 
probability and further enhance network connectivity [14]-[20]. 

Basically, the channel Rayleigh fading is conducted in systems introduced in the above-mentioned 
studies. In general, the dense obstacles and scattering are in the channels between the base station (BS) and 
the user which follow Rayleigh fading. Since the propagation encounters fewer obstacles and less scattering, 
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the Nakagami-m distribution is suitable to model the BS-relay links compared with Rayleigh fading. A key 
factor, channel fading makes a crucial influence on the performance of NOMA systems. In [21], the authors 
have studied the NOMA system over Nakagami-m fading to highlight outage performance, but the perfect 
channel estimation is assumed. Actually, the performance is further degraded by imperfect channel estimation. 
Therefore, it is very challenging to provide an accurate analysis of the NOMA system under joint impacts of 
imperfect channel state information (CSI) and Nakagami-m fading. 

Although the existing works about NOMA have assumed perfect CSI or the order of the instantaneous 
channel gain at the transmitter side [22]-[25]. However, this assumption is not suitable in some communication 
scenarios, such as underwater acoustic (UWA) [26] and high-speed railway (HSR) [27] systems, due to the 
rapidly changing channel and the large feedback delay. Motivated by such analysis, we consider NOMA 
system in imperfect CSI case. 


2. SYSTEM MODEL 

In this paper, we assume the base station(S), two user NOMA(U;,7 € {1,2}) shown in Figure 1. 
These nodes are equipped single antenna. To improve bandwidth usage efficiency, the user 1 (U,) works in 
full-duplex (FD) mode. Channels in link S-U1, S-U2, U1-U2 are denoted by h1, hz and hg, respectively. Self- 
interference channel existing at user U; is denoted by gı. These channels are followed Nakagami-m channel 
model. 





Signal link ——— > Interference link — — > 


Figure 1. System model 


The error channel estimation coefficients is modeled by [28] 
hj = hj +), (1) 


where h;(j € {1,2,3}) denotes the estimated channel coefficient and e; denotes the estimated channel error 
which can be approximated as a Gaussian random variable with CN’ (0, o?) 

In the block time ż, the base station S transmits the information axı (t) + a2%2(t) to U;. It is noted 
that xı and xə are the information of U; and U2, respectively To proceed NOMA, a, and az are the power 
allocation coefficients where a? + a3 = 1 and a? > a3. Thus, the received signal at U; can be computed by 


yu, (t) = VPs (azı (t) + azar (t)) (fu + e1 ) +4/Pigite (t=T) +n 


= 2 
= y Pshı (azı (t) + a2£2 (t)) + y Pseı (azı (t) + agX2 (t)) + VPigixv2 (t = T) + ni (2) 
ey 
effective noise 


where Ps and P; are the power transmission of S and Uj, respectively, n; is the additive white Gaussian noise 
(AWGN) with CN (0, No), gi is the residual interference coefficient and 7 denotes signal processing delay. 
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Then, the signal to interference plus noise ratio (SINR) to detect x2 at Uj is given by 
~ 2 
Ps a? | hy | 


Di2 = 





212 
Psa |ha| + Psoi + Pi Igal? + No 
a 3) 
mi 





x fd 
pa? [ña + po? + pig? +1 


where p = xe = 5- is the signal-to-noise ratio (SNR) and the SINR at U; to detect the own signal xı is given 


as 
alr 2 
ni 
r= — | (4) 
por +pl|m| +1 


Now, the base station S transmits signals to U2 through two links (direct and relay links). The received 
signal at U2 in direct link is expressed as 


yu, (t) = VPs (ax (t) + azz2 (t)) (he + e2) + ng 
= \/Pgho (azı (t) + azx2 (t)) Y Pse (a121 (t) + azz2 (t)) + n2 (5) 
eS SS 


effective noise 


It is noted that the signals received at Uz from U; node could be formulated as 


yu, >v (t) = /Pr (hs + es) zə (t — T) +n (6) 


where ng and ng are the AWGN with CN (0, No). 
Then, the SINR to detect x2 of U2 is given as 











fe 2 
oi 
Ty = ——, (7) 
paz ñ +poz+1 
From (6), the SINR to detect x2 at Uz is formulated by 
~ 42 
r P hsl 5 
2,1 — por + 1 . ( ) 
3. PERFORMANCE ANALYSIS 
3.1. Channel mode 
The probability density function (PDF) of their gains follow gamma distributions is formulated as 
= 9 
fine @) T (mj) (a) i K 
where m; and Q; are the fading severity factor and mean, respectively. 
3.2. Outage probability of U: 
The outage probability of VU, can be expressed as [29] 
OPu, = 1— Pr (T12 > h2, I1 > Yen1) (10) 
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where Ythi = 2Ri — 1 and R; is the target rate. With the help of (3) and (4), it can be written as 























~ 2 -2 
ti ni 
OPy, =1-Pr 3 ; > 2, > aa > 
paž [in| + poi + pial +1 poy + plat + ans 
_p ¢ (p0? + plal? + 1) 
=1-Pr ||| > 
p 
where ġ = max (we. 4). Then, (11) is calculated as 
1 2 2 
OPy, =1 - f har (x) i fp, aya (12) 
0 (po? +ex+1) 
Putting (9) into (12), we have 
1 m mar f mg Miki 
OPy, =1 ( 2.) Mia 
! T (Mg) T (m;,,) Qo, Qa 
* Sine ~S SN, (13) 
x feae Ro” f yte Wi“ dyda 
0 (po? +er+1) 
Based on [30, Eq. 3.351.2], (13) can be computed as 
m- (por mz; —1 k 
1 Mg me í NG i+) an 1 mg, (poz + px + 1) 
OPy, =1- ——~ e hy 
T (mg, ) Qg, k=0 k! Pia 
oo oe (14) 
P ety Pama A 
x feae “a1 "ha da 
0 
By applying [30, Eq. 1.111], (14) can be further simplified as 
mz; —1 bon om; (o2+1/p k eN 
TEN n a apoia ae ina Ar = to 
$ SN KIT (mg, ) Qz Qg 
oo (15) 





mg1%,, tom; 91 a 
x Jarci why dx 
0 


Finally, using [30, Eq. 3.351.3] the outage probability in closed-form for user U; can be obtained as 








em; (o24+1/p) mp, —1 REA n 7 7 
E a ) is 3 ( k ) r (Mg, +n) (o? + 1/p) (Qg) (omz, ) (Mg) õi 
l k=0 n=0 LM (7, ) k! (Op), E (mg, OQ, ae em, Ory ET 

(16) 


3.3. Outage probability of U> 
In this case, it is assumed that the far user U2 exploits the selection combine(SC) approach to process 
two paths achieved from the S and U1. Thus, the outage probability of U2 is given as [29] 


OPv, = Pr (max (T2, min (11,2, T'2,1)) < 72) 


=Pr (T2 < 2) Pr (min (T1,2,T21) < y2) (17) 
Ai Ao 
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Then, putting (7) into the first term of (17), we have 


Aj =1-—Pr(T2 > 2) 


2 
1p (iaf > BEE? ) 


p (af — 203) 

















oo (18) 
=1— i Fiial (x)dz 
y2(p02+1) 
p(aĵ = 245) 
Similarly, A, is obtained as 
ymy (eo3+1) m;,,—1 2 q 
Aya=1—e Paena So 1 f tome, (903 +1) a 
qo q! pQ;,, (aj — 7245) 
Furthermore, the second term of (17) is expressed as 
Ag = Pr (min (Ti2, T21) < y2) (20) 
=1— Pr (Ti,2 > y2) Pr (T21 > y2) 
Then, we can write Pr (T1,2 > 72) as 
2 R (po? + plo? + 1) 
Pr (T 1,2 > 42) = Pr | > 5) 5) 
(ay — V203) P 
oo oo (21) 
z= T fiol? (x) / Piat (y)dydz 
y2(eo?+el91|2+1) y2 (po? +pz+1) 
(a? -v2a2)P (a? -7v2a3)P 
Similar in above, Pr (T12 > y2) can be obtained as 
rami (ethe) mya g P (mg +n)! (o? + 1/0)” 
Peis gane Cy 5 2 (i la 7 Ta mg; 
=0 n=0 oi) k! (9.0%, ) (Mmg N7) (22) 
k -m <7} 
2M, ( VM; Qa, ) É 
x | ——— > Mor +- 
i — 7203) Ne aoh (af = 7243) 
Next, Pr (T2,1 > 72) can be formulated as 
oe Oe y2 poż +1) 
Pr (T21 > %2) = Pr (f 5 BI res fip (@ (23) 
Ta leegt) 
Similar, after some manipulations, we have 
a iy (eo3+1) Mi, —1 7 2 Pp 
a NNN} mz (pos +1 
Pr (Taa > x2) m PIi 5 Y2 hg (p 3 ) (24) 
are Pi 
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Putting (22), (23) and (26) into (20) and then with help the result and (19). Finally, the outage proba- 
bility of Uz is formulated by 





yams (p03 +1) Mj, —1 2 q 
OPa iae aey L (p03 +1) 
: q=0 q! PR, (a? E 7202) 








yom; (pod mz —1 
s hg cee Ymg, (p03 + 1) ” 
x [l—e E ye l 
= P! Pi 
(25) 


R 2 Š 
EA (c?+1/e) mi 1 k 


xe Ged Y > k ) I (mg, +o)" 
f= ee T (1mg,) k! (Qg N7, ) (Moa N7) a 


k = = 
x YM (m Q- + yom, Qg, ) Mg, =n 
(aj — 7243) Qh, SE m (aj — 7203) 


3.4. Throughput 
Based on achievable outage probability, we can compute throughput as below 


TP =(1-OPpy,)Ri + (1 — OPu,)Re (26) 





4. NUMBER RESULT 

In this section, we perform simulations to check the corrections of derived formulas. We set a? = 0.8, 
a2 = 0.2, Ry Rə 1, o? o? o3 o3 0.01, m = mg, mj, 
Q, = Q;, = 1 and Qo, = 0.05. 

Figure 2 shows outage performance of user U; versus the transmit SNR. We can see FD case exhibits 
its superiority compared with HD case. At high SNR region, SINR would be better and hence such outage 
probability can be improved at high SNR. Similar trend can be seen in Figure 3 for case of user U2. We change 
value of fading severity parameter m and a gap exists when comparing these cases of m. Similar performance 
is presented in Figure 4, Figure 5 for user U1, U2 respectively for two cases of noise term. Finally, throughput 
meets highest as SNR is greater than 20 dB shown in Figure 6. The reason is that throughput depends on outage 
probability. 
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Figure 2. The outage probability of U; versus SNR varying m with o? = 0.01. 
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Figure 3. The outage probability of Uz versus SNR 
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varying o? with m = 2. 
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Figure 5. The outage probability of U2 versus SNR 


5. 
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varying o^ with m = 2 


CONCLUSION 
This paper has studied cooperative NOMA scheme under impact of imperfect CSI. Based on the nature 
of cooperation, the near user benefits from full-duplex design to serve far user, which improve performance of 
far user. This results in an enhancement of the downlink NOMA over Nakagami-m fading. We analyze the 
system outage probability and throughput under different variations of the key system parameters. Finally, 
simulation results validate our considered system. 
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